Controlling coherence using the internal structure of hard n pulses 
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The tiny difference between hard -n pulses and their delta-function approximation can be exploited 
to control coherence. Variants on the magic echo that work despite a large spread in resonance offsets 
are demonstrated using the zeroth- and first-order average Hamiltonian terms, for 13 C NMR in Ceo- 
The 29 Si NMR linewidth of Silicon has been reduced by a factor of about 70, 000 using this approach, 
which also has potential applications in magnetic resonance microscopy and imaging of solids. 
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In magnetic resonance, a control pulse is hard if 
the pulse amplitude is much greater than the spectral 
linewidth and any resonance offset; hard pulses are often 
approximated as instantaneous delta- functions [li 13, [3| . 
The corrections to this picture are quite small for a single 
hard pulse, but they can lead to surprisingly large effects 
[J, |5| in important nuclear magnetic resonance (NMR) 
experiments that use many 7r pulses, such as the Carr- 
Purcell-Meiboom-Gill (CPMG) experiment @, 0j- Even 
though CPMG uses a very simple pattern of pulses [l(, 
Coherent Averaging Theory [8|] shows that the zeroth- 
and first-order correction terms arising from non-zero 
pulse duration are quite complicated [4J, [5(, making a 
quantitative prediction of their effects very difficult. 

In this Letter, building upon our earlier results [4J, l5(, 
we design more complicated pulse sequences, and show 
that much simpler approximate Hamiltonians can quan- 
titatively explain the experiments. This shows that the 
small difference between hard it pulses and their delta- 
function approximation can be put to good use, enabling 
new classes of spin echoes which have promising appli- 
cations in NMR, magnetic resonance imaging (MRI) or 
microscopy of solids, and related spectroscopies. 

The NMR data in Figs. 1-4 of this paper were ob- 
tained with powder samples (Cgo or Silicon doped with 
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Sb (10 /cm 3 )) at room temperature, in B ext =12 Tesla. 
Both samples are well-approximated as a single species 
of spin 1=1/2 nuclei ( 13 C or 29 Si), coupled together by 
the like-spin dipolar interaction [J, |5|. For a meso- 
scopic cluster of N-spins, the Hamiltonian in the rotating 
frame is Hint = Tiz + 7~L Z z , where a net resonance off- 
set ( Sl™ et =Vl 9 ?r° iet + Q l ° c \ gives rise to the Zeeman term 
TLz=Q z let I ZT1 and the secular part of the homonuclear 



J,,-h-h 



dipolar coupling P,0,[3 is H Z! ^2f >i Bij(3I ZiJ . Zj -^-^ ; 
Our macroscopic powders are similar to an ensemble of 
N-spin clusters, with distinct fi z ° c values in different clus- 
ters due to bulk diamagnetism |5j . The resulting Zeeman 
line broadening dominates the spectrum's full width at 
half maximum (FWHM), which was only about 2 ppm 
(e.g., the 13 C ( 29 Si) spectrum's measured FWHM = 260 
Hz (200 Hz), while the calculated dipolar FWHM=38 Hz 
(88 Hz)). The rf pulses used were unusually hard (e.g., 



FIG. 1: (color online) Sample Ceo- (a) Comparison of 
CPMG (grey) to APCPMG (blue). Inserting a single flip- 
180y pulse into APCPMG induces an echo of the echo 
train (green), (b) Inserting a single 180x (red) has no ef- 
fect (blue). (c) Reversing the APCPMG phase pattern, 
90x - {-Y,YY°° ~{Y,-Y} 600 , at the point indicated has 
the same effect (black) as inserting a single 180y pulse 
(green), (d) A CPMG of the echo train is induced by using 
90x~{-Y,Y} w -({Y,-Y} 20 -{-Y,Y} 20 ) repeat . For (a-d), 
t = 25fis, f2|C°»" et = 0, a ~ 0.71, and only the peak of each 
echo is shown. The signals in Figs. 1-4 are normalized to the 
amplitude of the Ceo or Si:Sb FID signal. 



the pulse strength uji/2tt w 25 kHz (16.4 kHz) was about 
100 (82) times the 13 C ( 29 Si) linewidth, with a 128.56 
MHz (101.56 MHz) Larmor frequency 4]). Low coil fill- 
ing factors [1] (<8% for 13 C data and -40% for 29 Si data) 
made the rf pulses very uniform across the samples. 

The open grey squares in Fig. QJa) show the amplitude 
of each peak in a long-lived train of spin echoes [9( gen- 
erated by the Carr-Purcell-Meiboom-Gill (CPMG) ex- 
periment, 90^ — {Y, Y} , where the first pulse is a 90° 
rotation about the X-axis in the rotating frame [y, |7|. 
The block {Y, Y}, repeated N times, represents the se- 
quence (r - 180y -It - 180y - r) where the 180° ro- 
tations are about the Y-axis, and echoes are acquired 
in the 2r time interval after every 180° (or it) pulse 
[li \A \M- I n contrast, the train of spin echoes quickly 
decays to zero (Fig. HJa), blue) for Alternating-Phase 



CPMG (APCPMG) 90 x - {-Y,Y} N . 

To understand this dramatic difference, we apply Co- 
herent Averaging Theory [8] to the repeating block 
{</>!, </>2j-, with f80° pulses of duration t p about the </>i 
or cf>2 axis, and cycle time t c — At + 2t p . Short t c is 
used throughout this paper, so it is a good approxima- 
tion to keep just the first two terms H^ + TiS 1 ' in the 
Magnus expansion [|. The {Y, Y] block has [| H^ F} = 
afi zz — f3H y y=Il, while the {— Y, Y} block has a slightly 



different form: 
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for <t=x, y, or z. The extra term — AH,™ et I XT looks like a 
constant transverse field in the X-direction, which, when 
acting alone, causes spins to nutate 1] in the Y-Z plane in 
a manner we define as clockwise (CW). Variation in £T" et 
values across the macroscopic sample leads to a spread 
in precession angles that causes signal decay. In the 
well known free induction decay (FID), T% arises from 
a spread in fl" et of the original Zccman Hamiltonian. By 
analogy, the rapid decay of the spin echoes produced by 
90 x - {-Y,Y} N (Fig. []Ta), blue) can be thought of as 
an 'FID of the echo train'. 

Attempting to undo this XJ-like decay, we insert a 
single 180y pulse into the APCPMG sequence, 90x — 
{-Y, Y} Nl - 180y - {-Y, Y} N \ which produces a strik- 
ing 'echo of the echo train' (Fig. HJa), green). Al- 
though this looks like a conventional Hahn echo [9(, the 
signal actually extends over more than 800 individual 
spin echo peaks. The dephasing caused by —\Q™ et I XT 
(CW precession) during the Aq£ c is followed by counter- 
clockwise (CCW) precession caused by +\fl™ et I XT , and 
this rephasing leads to the echo of the echo train when 
Nz=Ni. When a single flip-180x is used instead, no echo 
of the echo train (Fig. QJb) red) is seen, as predicted by 
our model, because a perfect rotation along the X-axis 
does not change the sign of the — A£l™ et I XT term. On 
the other hand, the echo of the echo train (Fig. QTc) 
black) is recovered if the flip-180x is removed and the 
phase pattern in the second repeating block is reversed 
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from {-Y, Y}to {Y, -Y}, since [5|] H^_ y} =H&0%-i x t, 

compared to 1if > _ YY} =E.-AQ,™ t I xT , In Fig. [^c), the 
phase reversal of 1200 hard 7r pulses yields a signal indis- 
tinguishable from that induced by the single fhp-180y, as 
predicted by our model. In contrast to this model, taking 
the limit of delta-function pulses (t p — > ) would kill [J,l5| 
the transverse field terms in H^+H^ exploited here and 
throughout the rest of the paper. FigureQJd) shows that 
the approach of Fig. |Tfc) can be repeated, creating mul- 
tiple echoes in the envelope of individual spin echo peaks, 
or a 'CPMG of the echo train'. However, the signal does 
decay, since the sign of the term H is never reversed in 
Fig. HJ To beat this decay, we use an approach inspired 
by the magic echo [l|, [lfj, 111 | • 

In the original magic echo 1CJ, [ll|, [l2| , a continuous rf 




FIG. 2: (color online). Sample Ceo- Three experiments 
inspired by the magic echo [10L [tH |. which all start with 
{— X, X} , have distinctly different results. With 90-x fol- 
lowing the repeating block, no magic echo forms (red); with 
90+x following the repeating block, a large echo emerges 
(blue) ; when applying a 180y pulse at time tf ± after the burst 
of the failed sequence (red), an optimized echo is achieved 
(green). Here, N=200, r = 50^s, 0^°^= 0, and a w 0.83. 



field in the transverse plane picks out the part of the dipo- 
lar coupling that is secular in the strong transverse field 



0. In the {-X, X} block d, the effective field AO™ e *I, 



in H 



(0) 

{-x,x\ 



VT 



could play the same role, as first proposed by 



Pines and Waugh 13] for a single value of Sl™ e *. Figure 
[2] provides experimental support for their prediction, 
even though the weakness of the effective transverse field 
makes it hard to justify the second averaging analysis 
[la . Il4 |. In addition, the spread in $7™ e * across the 
macroscopic sample has non-trivial consequences, as 
shown by the different effects (Fig. [2]) of the two bursts, 
{— X, X} N — 90±j>f, followed by a free evolution of du- 
ration tf re e- Using our model, the unitary operators are 

e -i(W„+nr t ^T)*/Vee e -i(^^W„TAnr t /,T)A r tc Wgo±xj 

where ±X is the 90° pulse phase and a > (3 for our 
experiments [151 ] - 

For the -X choice, the Zeeman phase wraps in a CCW 
manner both during and after the burst, which spoils 
the magic echo that would otherwise form during the 
free evolution period (Fig. red). For the +X choice, 
both Zeeman and dipolar terms switch from CW phase 
wrapping in the burst to CCW phase unwrapping dur- 
ing the free evolution period, resulting in a large echo 
(Fig. [2l blue). This echo is not optimized, since the 
rcfocusing time is different for the dipolar and Zccman 
phases {t dipo i ar = (a - 0) Nt c /2, tzeeman = XNt c ). An 
optimized echo (Fig. [21 green) is generated if we apply 

a 180y at time t fl = ( °~^~ 2A ) Nt c after the failed se- 
quence {-X, X} N - 90_x (Fig. H red). This sequence 
aims to synchronize the refocusing times of the dipolar 
and Zeeman phases by using the fact that a 180y pulse 
flips the sign of the Zeeman term but does not change the 




To test our model, we use phase-coherent frequency 



FIG. 3: 



(color online) . 



Sample Ceo- 



Using 



fl l ° c ± ^offset' the quadratic echoes produced 

JV JV 

by {X,X}~{-X, -X}- - 90y - t free : green (^offset = 
Hz), black (v ff se t = —3 kHz), differ from the linear 

N N 

echoes produced by {— X, X}~2 {X, — X}~2~ — 90x — tfree- 
blue {^offset = Hz), red {u ff set — -1 kHz), where 
^offset ~ ~hv ffset- Only the black echo shifts to the 
right. (Inset) Image plot of 31 quadratic echoes for Hz 
< 0. 3 of ob f f e Jh < 3 kHz, in steps of 100 Hz. The black trend 
line shows our predicted Zeeman refocusing time. Here, 
N=I00, r = 10/US, and a « 0.5. 



dipolar term. The measured echo happens at a slightly 
different time, due to terms ignored in this model 15]. 

Compared to the original magic echo [lfj, [ll|, which 
works best if fi™ e *=0, sequences based on the {— X, X} N 
block have several clear differences: both Zeeman and 
dipolar phases are wrapped during the burst, a 90±x is 
used instead of the 90y , and the 2r gaps in between the 
■k pulses of the {— X, X} block simplify implementation. 

The {</>i, </>2 } blocks used so far have an effective 



transverse field term in TL 
>(0) 



(0) 

{01,02} 



However, for the block 



{X, X}, TL, X x ^=aTL zz —/3TL xx , and so the first transverse 



L {x,x} 
field term is in TL 



(i) 

{X,X} 



yn\t, z ) 1 x T ' ri -^x.X} 



where k 2 = t p (8r + 2t p ) / (2t c fnr) 5]. In principle, 
despite its smaller size, the quadratic transverse 



field term of TL 



(i) 

{X,X} 



could be exploited just like 



(0) 

{-x,xy 



the linear transverse field term found in TL 
In practice, however, {X, X} N is a poor nutation 
experiment, since ft, x ' x , XT causes rapid signal 
decay. Inspired by the Rotary Echo experiment [16l |. 



,N 



we tried replacing {X, X} with the composite block 
{X,X}*{-X,-X}^" because Hf]_ x _ x} = TL {0) 



and TL 



(i) 

{-X-X} 



{X,X} 



-7~Lr x x\ H, an d we managed to 
recover most of the original signal. We thus infer 

I L N N 

EH that the net effect of {X,X}^{-X, -X}^ is 
well-approximated by the much simpler unitary operator 

where we allow for different Q™ e ' during {±X, ±X} 2 . 



■>nei,± . 



-n l r ± ^fYL (WjftLm during the 



jumping^ -,, z _l. " offset v "off s 

burst {X, Xp{-X, -X} f -9(V, followed by n 9 o l f o *f et =0 
during free evolution, leading to [15(] the unitary operator 



■^global 



oc h T )Nt cJ , 



Increasing £l ri set increases the Zeeman dephasing dur- 
ing the burst, pushing the quadratic echo peak out 
to later in tf ree (Fig. [3l green and black). The inset 
of Fig. [3] shows the strong agreement between the 
Zeeman refocusing time predicted by our model (black 
trend line) and the quadratic echo peak measured in our 
experiments over a range of ^ 9 t°f^ et - In contrast, the cor- 

JV JV 

responding linear sequence {—X, X} 2 {X, —X} 2 — 90x 
has its largest signal just after the burst, for all Sl 9 „?t" et 
(Fig. [31 blue and red), as predicted in our model [15| . 

Controlling both dipolar and Zeeman phase wrapping 
using TL, 1 . i is an unusual aspect of the quadratic 
echo. As one use of this, we designed a composite 
block with no net dipolar evolution over duration of 6A, 

JV JV 

(A-h5)-90 V)1 -{X,X} T {-X, -X} T -90 fe -(A-<5), which 
we refer to as {N,S, ^1,^2}, with A = Nt c /4, \S\ < A, 
and ipi — i^" f° r i=l,2. For constant fi™ et , the uni- 
tary operator is U\^ Y e h[ - * "t)^~<~ > for tjj 1 = ip 2 , and 
e -|(nj e */ ZT )(+2A) for ^ ^ ^ rj]| While similar effec _ 

tive operators were previously demonstrated [17j | using 
magic sandwich echoes for ||Wz|| < C]|'Hza||, our approach 
works in the complimentary regime ||Hz||>||H Z z||, where 
the scales are calculated using [2J, [5j || A|| 2 = Tr(A^A). 
In particular, the {X, S, ^1, ^2} sequence is still effective 
even when there is a large spread in £7™ e * values across 
the sample [15j. 

Our model predicts that both Zeeman and dipolar 
phase are refocused after each fJVO, tpi, tpi} block, yield- 
ing a time-suspension sequence [l7(. Indeed, in Si:Sb, our 
sequence pushes the decay time from T£ w 1.6 ms out 
to T2 c lve ~ 110 seconds, or about 10 10 periods of 
Larmor precession (Fig. 0Ja), blue), quite close to the 
spin-lattice relaxation time, T\ — 290 seconds. The nor- 
mal linewidth is thus reduced by a factor of about 70,000 
((Fig. Ha), Inset). 

Eliminating dipolar dephasing in order to measure 
O" 6 * in applied magnetic field gradients enables the MRI 
US El, EH or MR microscopy (2(j of solids. Measuring 
the spectrum in a field gradient is the first step toward 
imaging using the back-projection technique [lj.ll9|. Fig- 
ure [UJb) shows a faithful reproduction of an input top- 
hat spectrum, where each spectrum is the Fourier trans- 
formation of the pseudo-FID resulting from two inter- 



laced data sets [15j. Note that both the signal ampli- 



tude and the v ffset values have been quite accurately 
reconstructed in this approach. Compared to existing 
approaches for the MRI of solids 17], |l8|, Il9| . our ap- 
proach does not need to switch off the applied Zeeman 
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FIG. 4: (color online), (a) Sample Si:Sb. The 29 Si time- 
suspension data using the sequence 90x — {2, 0, —Y, — y} 84000 
with t = 60^is, v ffset~ 2.5 kHz (blue) and corresponding 
fitting curve (black), extend far beyond the normal 29 Si FID 
with v ffset= Hz (red), (a, Inset) The 200 Hz normal spec- 
trum (red) is narrowed to 0.003 Hz (black, Fourier transfor- 
mation of the fitting curve), centered at v ffset- (b) Sam- 
ple Ceo- Reproduction of a top-hat lineshape using sequence 
90x - {2, to, -y, -Y} - {2, 0, -y y} 30 with r = 22^s, and 
io=0. Each trace is the measured spectrum of a pseudo-FID 
with different v ffset , f° r -4 kHz< v ffset <+4 kHz in steps 
of 500 Hz, covering the range 2ir|f //»et[/wi < 16%. To ob- 
tain this full bandwidth, the pseudo-FID interleaves a second 
data set using the same sequence, but with to — — (-y + r' - )• 



gradient inside the bursts, which enables the application 
of large field gradients at moderate cost. It should also 
be possible to implement standard frequency- and phase- 
encoding methods using this approach [l|,ll9(. Since pulse 
strength varies across a big sample, the uniform pulse 
assumption of our model is a potential concern. Exper- 
imentally, an intentional uniform misadjustment of all 
pulse angles leads to similar MRI top-hat lineshapes and 
to similar line-narrowing performance, suggesting that 
these two sequences are robust 15]. 

Our sequences may help in the study of some impor- 
tant biomaterials, since the Hint assumed here is very 
similar to that of 31 P in bones and teeth 
inary results are encouraging 
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Prclim- 
15j. These sequences also 



have potential applications in proton ( 1 H) NMR. While 
the dipolar linebreadth dominates most H spectra, a 

large ^offset can ^ e use< ^ to reacn the || Hz \\>\\ H zz \\ 
limit of our model, as demonstrated in our preliminary 
results on Adamantane [151 ] . Future work will use micro- 
coils [23|, |24J, |25j to reach shorter t c , which should improve 
the utility of our model for proton NMR experiments. 

Related effects can occur for a wider variety of Hint 
and Hp^ than we have treated here, provided that 
[Hp^ , Hint] 7^ 0. Shaped pulses, soft pulses, and strongly- 
modulated pulses have proven to be important elements 
of the NMR toolbox. Exploiting the internal structure of 



hard ir pulses provides us with yet another technique to 
control the coherent evolution of quantum systems. 
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